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Abstract. The model of inelastic Maxwell particles (IMP) allows one to derive some exact results which show the strong 
influence of inelasticity on the nonequilibrium properties of a granular gas. The aim of this work is to propose a simple model 
kinetic equation that preserves the most relevant properties of the Boltzmann equation (BE) for IMP and reduces to the BGK 
kinetic model in the elastic limit. In the proposed kinetic model the collision operator is replaced by a relaxation-time term 
toward a reference Maxwellian distribution plus a term representing the action of a friction force. It contains three parameters 
(the relaxation rate, the effective temperature of the reference Maxwellian, and the friction coefficient) which are determined 
by imposing consistency with basic exact properties of the BE for IMP. As a consequence, the kinetic model reproduces the 
true shear viscosity and predicts accurate expressions for the transport coefficients associated with the heat flux. The model can 
be exactly solved for the homogeneous cooling state, the solution exhibiting an algebraic high-energy tail with an exponent in 
fair agreement with the correct one. 
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INTRODUCTION 

The prototype model for the description of granular media in the regime of rapid flow consists of an assembly of 
(smooth) inelastic hard spheres (IHS) with a constant coefficient of normal restitution a < 1. In the low density 
limit the velocity distribution function obeys the Boltzmann equation (BE), modified to account for the inelasticity of 
collisions 1 1 ]. Because of the mathematical intricacy of the BE for IHS, a simpler model of inelastic Maxwell particles 
(IMP) has been proposed Q> S- Si - where the collision rate is assumed to be independent of the relative speed of the 
colliding pair. Apart from its interest as a model of granular gases, the IMP model is interesting by itself since it 
allows the derivation of some exact results for any dimensionality d, showing unambiguously the strong influence of 
inelasticity on the nonequilibrium properties of the gas. 

While the BE for IMP is considerably simpler than for IHS, it is still a formidable task to solve it in a closed form, 
even in the case of the homogeneous cooling state (HCS). The aim of this work is to propose a simple model kinetic 
equation that preserves the most relevant properties of the BE for IMP and reduces to the celebrated Bhatnagar-Gross- 
Krook (BGK) model in the elastic limit a — > 1 . In the proposed kinetic model the collision operator is replaced by a 
relaxation-time term toward a reference Maxwellian distribution plus a term representing the action of a dissipative 
friction force. The kinetic model contains three parameters: a relaxation rate modified by a factor /3(a) with respect 
to its elastic value, an effective reference temperature modified by a factor 9(a) with respect to the actual granular 
temperature, and a friction coefficient y(oc). The model is a hybrid of two previous models |5, 6] originally proposed 
for IHS , reducing to them if either y( a ) = or 9 ( a ) = 1 , respectively. The three parameters J3 (a), 9(a), and y( a ) are 
determined by imposing consistency with basic exact properties derived from the BE for IMP. As a consequence, the 
model reproduces the true shear viscosity and predicts accurate expressions for the transport coefficients associated 
with the heat flux. Moreover, it can be exactly solved for the HCS, the solution exhibiting an algebraic high-energy 
tail with an exponent in fair agreement with the correct one, especially for high inelasticity. 



THE BOLTZMANN EQUATION FOR INELASTIC MAXWELL PARTICLES 

The BE for IMP can be obtained from the BE for IHS by replacing the term |g- a\ in the collision rate (where 

g = V] — \2 is the relative velocity and a is the unit vector directed along the centers of the two colliding spheres) by 
an average value proportional to the thermal speed \J2T jm (where T is the granular temperature and m is the mass 



of a particle). In the version of the IMP model first proposed by Bobylev et al. |2] the collision rate has the same 
dependence on the scalar product g • G as in the case of hard spheres. In a simpler version 01 0L the collision rate is 
assumed to be independent of g • a. In this latter case, the corresponding BE reads 

V)/(r,T i; =J[w,t\f\ = d t 2 V f V :l [da[ dv 2 (a-'b- 1 - l)/(r,v i;f )/(r,v 2 ;f), (D 

2 n(r,t)Ll d J J 

where n is the number density, Vo « nT x l 2 is an effective collision frequency whose specific form will not be needed, 
Cl c i = 2n d l 2 /T(d /2) is the total solid angle in d dimensions, and b is the operator transforming pre-collision velocities 
into post-collision ones, namely, bvi, 2 = Vi, 2 T (1 + °0(g ' 0)a/2. The Boltzmann collision operator J[f] conserves 
mass and momentum but not energy. The collisional moments of second and third degree are |7| 

'^JdyV 2 J[f]=-^p, mJdy(V i Vj-d- 1 V 2 8 i j)j[f} = -Vr,(P i j-p8 i j), ^ J dyV 2 VJ{f] = V K q. (2) 

Here, V = v — u is the peculiar velocity, where u is the flow velocity, P is the pressure tensor, p = nT = d~~ trP is the 
hydrostatic pressure, and q is the heat flux. The exact expressions for the cooling rate £ and the rates of change v-q and 
v K are 01 



(4) 



In isotropic states, the collisional moment of fourth degree is 

n- 1 J dvV 4 J[f] = -v 2 M 2 + Av (2r/m) 2 , M t = rT l J d\V 2t f, 

where the exact expressions for the rate of change V2 and the coefficient X are |0] 

v 2 * = v 2 /v = {l + a) 2 (4d-l + 6a-3a 2 )/\6d + 2^* 1 X = (1 + a) 2 (d + 2)(4d- 1 -6a + 3a 2 )/64. (5) 

In the case of the uniform, free cooling state, the BE Q and the evolution equations for the second- and fourth- 
degree moments become 

d,f{v)=J[v\f], d t T = -£T, d t M 2 = -v 2 M 2 + Xvo(2T/m) 2 . (6) 

The solution to the cooling equation is T(t) = r(0)[l + ^(0)f/2] _2 (Haff's law). If time is measured by the ac- 
cumulated number of collisions per particle as x{t) = f Q dt' v Q {t') = (2/£*)ln[l + C(0)f/2], Haff's law becomes 
T(z) = r(0)exp(— It is convenient to introduce the reduced moments M\ = Mi/(2T/m) e , so that M\ =d/2 
by definition, and the reduced distribution f* (c, t) defined by 

f{v,t)=n[m/2T{t)} d > 2 f{c,x), c = v/y/2T(t)/m. (7) 

Thus, Eq. l|6) yields 

d*f{c) + (C/2)d e .cf(c)=J*[c\f], d x M* 2 = -{v* 2 -2t;*)M* 2 +X, (8) 

where 7*[c|/*] is the dimensionless version of the collision operator /[v|/]. Except in the one-dimensional case, the 
reduced moment M\ (t) converges in time to the well defined value M\ (°°) = A/(v|— 2f *). In general, the distribution 
function reaches a scaling form 0, called homogeneous cooling state (HCS), which is the stationary solution of 
Eq. (|8}, i.e., /*(c, t) — > /*(c,°°) = /hcs( c )- ^ s exact expression is not known, except in the one-dimensional case, 
where / h * cs (c) = (2 3 / 2 /tt)(1 +2c 2 )~ 2 0. For d > 2, the fourth cumulant (orkurtosis) a 2 = 4M^/c/(c/ + 2) - 1 of the 
distribution /J cs is 

^ 4 X 1= 6 ( 1 -«) 2 (9) 

2 d(d + 2)v^-2^* 4d-l + 6a~3a 2 ' 

Therefore, a 2 > for IMP, in contrast to what happens in the case of IHS with a > \[2j2 flflll . It is also known that 
/hcs ex hibits an algebraic high-energy tail of the form 

ti cs (c)~ c - d -« a \ (10) 



where the exponent s(a) is the solution of the transcendental equation 

, l ~ a \ ./o n ^2/,i , fl + a\ s r(s/2 + l/2)r(d/2) 



2^1 [a,fe;c;z] being a hypergeometric function. Equation dlOi implies that those moments M| with £ > s(a)/2 are 
divergent. 

Finally, the exact expressions for the transport coefficients are |7] 

_ p 1 _d + 2 p l+2a2 _2T , ,_d + 2 p l+3fl2/2 

r,_ VbV*-C*/2' K_ 2 mv V*-2C*' M_ T i,C_ J ' * ~ 2 mv V*-3C*/2' U J 

Here, 77 is the shear viscosity, K is the thermal conductivity, k' is a modified thermal conductivity, and pi is a coefficient 
relating q and Vn. As can be seen from vjj > £*/2 and v£ > 3£*/2, so that 77 and k' are well defined for all a. 
On the other hand, K and ^ are not positive definite if v* < 2£*, i.e., if a < a/, = (4 — d) /3d. Therefore, if d < 4 and 
0! < Ct/,, there is no hydrodynamic behavior since the heat flux does not relax to a Fourier-law form. 



MODEL KINETIC EQUATION 

Although the BE for IMP is more manageable than for IHS and some important properties are accessible in an exact 
way, its explicit solution is not known, even for the HCS. It is then natural to wonder whether a simple generalization 
of the well-known BGK model kinetic equation to the case of IMP can be proposed. The model considered in this 
paper is constructed by the replacement 

y[r,v;f|/]^7[r,v;f|/] = -j3(a)vo(r,f)[/(r,v;f)-/o(r,v;f)] + 7 (a)vo(r,r)5vV/(r,v;f), (13) 

where 

/o(r,v;f) = n(r,t) [m/2ne(a)T(r,t)} d/2 exp [-mV 2 /29(a)T(r,t)] (14) 

is a local equilibrium distribution parameterized by the temperature 6(a)T. The effect of the inelastic collisions in the 
original BE is played in the model dl 31 by a relaxation term toward the distribution /o at an effective temperature 6T, 
plus an external friction term. The model contains three free parameters: the factor /3 (a) > modifying the relaxation 
rate with respect to its elastic value, the factor 9(a) < 1 modifying the granular temperature in the reference state 
/o, and the friction coefficient 7(a) > 0. These three parameters will be chosen in the next section to optimize the 
agreement with the most important properties of the BE for IMP. If one particularizes to 7(a) = 0, the model d 1 3I > 
reduces to the one proposed by Brey, Moreno, and Dufty (BMD) |5], while the choice 0(a) = 1 yields a simplified 
version of the model proposed by Brey, Dufty, and Santos (BDS) |6]. From that point of view, the model d 1 3I > can be 
seen as a hybrid of the BMD and BDS models. Although the two latter models were originally proposed for IHS, Eq. 
dl 3i is proposed here as a model for IMP, not IHS. In fact, it has been shown fTllfl2Tl that the BMD model shares more 
features with the BE for IMP than for IHS. 

Let us now obtain the basic physical properties of the model operator J. First, it is straightforward to check that 
it conserves mass and momentum. However, since 0^1 and 7 7^ 0, energy is not conserved by collisions. More 
specifically, the cooling rate and the rates of change defined by Eq. (0 are given (in reduced units) by 

C*(o)=P(o)[l- 0(o)] +27(0), v* (a) = ,3 (a)0 («) + £*(«), v* K (a) = ^[39(a)-l] + 3 -C(a), (15) 

where henceforth a tilde means that the corresponding quantity has been evaluated with the model operator J. 
Moreover, 

«-v/*v«^ =-(^w+/jM) (^y. (.6, 

In particular, setting I = 2 we reobtain Eq. @ with 

vf(a) = j3(a)[20(a)-l]+2£*(a), 1(a) = P(a)9 2 (a)d(d + 2)/4. (17) 



Let us now consider the free cooling problem. As shown by Eq. ©, the necessary and sufficient condition to reach 
a scaling solution (HCS) with a finite fourth-degree velocity moment is v| > 2£*, i.e., > 1/2. In that case, the first 
equality of Eq. l|9} implies that the fourth cumulant of the scaling solution is 

52(a) = [1- 0(a)] 2 / [20(a) - 1] . (18) 

Therefore, we have > 0, regardless of the precise values of the parameters /3, 0, and 7 of the model. As a 
consequence, the model dl 31 can never capture the negative values exhibited by in the case of IHS for a > V2/2. 
From Eq. (I16> it is easy to obtain the time-dependence of the reduced velocity moments in the free cooling problem: 

MKr)=e-^'-^^ (19) 

Since, on physical grounds, < 1, it turns out that the reduced moments M|(f) diverge in time if I > (1 — 0) _1 . This 
implies that the scaling solution presents a high-energy tail of the form JlQi with the exponent 

7(a) = 2[l-0(a)] _1 . (20) 

Finally, the transport coefficients of the model are given by d!2i . except for the obvious replacements — > £*, 
V* -> V*, V* -> V*, and a 2 -> a 2 . 

The main advantage of a kinetic model lies in the possibility of obtaining the explicit form of the velocity distribution 
function. Let us illustrate this in the free cooling case. According to the model Jl 3i . the first equation of l|8} becomes 

{ i 8- 1 «9 T + l + [(l-0)/2]«9 c -c}r(c)-(^0)- rf / 2 exp(- c 2 /0)- (21) 

It is interesting to note that the parameter 7 does not intervene explicitly in Eq. d2TY so that it is formally equivalent to 
the equation obtained from the BMD model [ 5] . Given an arbitrary initial condition /* (c , 0) , the exact solution to Eq. 

(ED is OB 

f*{c 1 T)=e-^ +d ^W\r (e-^ 1 - e ) T / 2 c,o)+(7T0)- c/ / 2 J PT dyex P {-[l+d(l-9)/2}y-c 2 e-^y/e}. 

(22) 

The HCS is obtained taking the limit T — > °° with the result 

/ h * cs (c) = (7r0)-'/ 2 (l - 6)-\e/c 2 ) d l 2+l ^-V [r (5/2+ 1/(1 - 0))-T (5/2+ 1/(1 - 0),c 2 /0)] , (23) 

where the change of variable x = c 2 exp[— (1 — d)y]/0 has been made and T(z,a) = /J° dxx z ^ 1 e^ x is the incomplete 
gamma function. Note that the whole dependence of yjL on inelasticity appears through the parameter only. In 
the high-energy limit (actually, if c 2 > 0), T (5/2 + 1/(1 - 0),c 2 /0) — > 0, so that we recover the tail (O with the 
exponent ( I20L The transient from the initial distribution /* (c, 0) to the asymptotic distribution /^ cs (c) is described by 
Eq. (I22> . which can be rewritten in a simpler form by introducing the deviation 5/*(c,t) = /*(c,t) — /h cs (c). The 
results for 8f* and its moments 8M* t with I < (1 — 0)~' are 

8f(c,T) = e-W + < 1 - e y*8f (e-W-'WcO) , 5M|(t) = e-^-^-^S^O). (24) 



While /*(c,t) relaxes uniformly to f^ cs (c) at fixed c with a relaxation rate j8[l +5(1 — 9)/2], the moment M|(t) 
relaxes to its HCS value with a shorter relaxation rate j8 [1 — 1(\ — 6)] which goes to zero as I approaches the threshold 
value (1 — 0)~' from below. This paradoxical property flUl is a consequence of the non-commutability of the limits 
c — » 00 and t — > °° in Eq. J24L To analyze this with more detail, let us suppose that the initial distribution /*(c,0) 
has a high-energy tail much weaker than that of yj cs (c), so that 8f* (c,0) ~ — /h cs (c) ~ — Ac~ d ~ 2 /( l ~ 9 " > in the region 
c 2 3> 0, where the amplitude A is known from Eq. i23\ . In that case, at fixed x, one has 



8f(c,T)^-Ac- d - 2 ^- d ^-f^(c), c 2 »0 e ^ 1 - e ) T . (25) 

This means that, at any fixed time T, there always exists an infinite range of large speeds, c 2 3> 0e^''~ e ' T , where the 
deviation of the velocity distribution function from its asymptotic HCS value is of 100%. This region is eventually 
responsible for the divergence of moments with £ > (1 — 0)~'. 




RESULTS AND COMPARISON WITH THE BOLTZMANN EQUATION 

So far, all the properties of the kinetic model (1 1 31 described in the preceding section are valid with independence of the 
specific expressions for the parameters j3 (a), 9(a), and y(a). Now we fix them by requiring the model to reproduce 
the basic properties of the original BE. The most characteristic consequence of inelasticity is the cooling rate, so that an 
obvious requirement is £* = £ *. The next requirement could be the agreement with either the relaxation rate v^j or v£. 
We cannot enforce both since that would imply 9 = d/(d + 2), which yields an unphysical model in the elastic limit. 
This impossibility of satisfying the shear viscosity and the thermal conductivity simultaneously, which also happens 
with the conventional BGK model, is remedied by a more sophisticated model 11211 . In the case of the kinetic model 
Jl 3i . let us take v^j = as a second condition. Finally, given the important role played by the kurtosis a 2 of the HCS, 
the third condition adopted here is a 2 = These three requirements yield 

^ (a) = ^p 7(«) = ^(l-« 2 )-^^[0-'(«)-l], 9{a) = l+a 2 {a)-^a 2 (a)[l+a 2 (a% 

(26) 

where a 2 (a) is given by Eq. (|9j. 

By construction, the model reproduces the exact shear viscosity r\ but not the transport coefficients associated with 
the heat flux. However, as Fig.n s h ows i the model captures reasonably well the rapid increase of fc(oc), K'(a), and 
H(cz) with the inelasticity, the agreement being especially remarkable in the case of ju(oc). The model also predicts 
the existence of a threshold value a/, below which hydrodynamics no longer holds. The value of a/, is the solution of 
a quartic equation which for d = 3 yields a/, ~ 0.1 14, in excellent agreement with the exact result a/, = 1/9 ~ 0.1 1 1. 

As Fig.[2]shows, the model underestimates the exact exponent s(a) but the agreement improves as the inelasticity 
increases. The time evolution of the distribution function in the free cooling state, as well as the asymptotic scaling 
solution are also shown in Fig.|5]for a few representative cases. 



CONCLUDING REMARKS 



The model of IMP shares with the more realistic model of IHS the description of the important influence of inelasticity 
on the dynamical properties of a granular gas. However, that influence is magnified by the IMP model, giving rise 
to stronger departures from the Maxwellian distribution and even to the absence of hydrodynamics for sufficiently 
inelastic systems. Nonetheless, the fact that non-trivial exact results can be derived from the BE for IMP justifies its 
study in order to gain a broader perspective on the peculiar properties of dissipative gases. From that point of view, 
the model Jl 3I > proposed here can be useful to have access, at least at a semi-quantitative way, to relevant information 
(such as the velocity distribution function itself) not directly available from the BE for IMP. 

The key ingredient of the kinetic model dl3t . also present in the BMD model |5], is the effective temperature 
9(a)T < T in the reference distribution fy. Its existence leads to a non-trivial HCS solution with a positive definite 
kurtosis and an algebraic high-energy tail. The additional presence of the friction term relieves the effective tempera- 
ture of fully accounting for the cooling rate, so that a closer contact with the BE for IMP is possible. 

The analysis presented in this paper can be extended along a number of routes. The explicit solution of the kinetic 
model to other problems, such as the simple shear flow, is straightforward. Moreover, the flexibility of the model 
allows one to choose its parameters by a fit to other quantities, such as the effective collision frequency V K (a) and/or 
the exponent s(a), different from the ones considered in this paper. The Gaussian kinetic model proposed in Ref. fl^l . 
which is able to reproduce V-q (a) and V K (a) simultaneously, can be extended by the inclusion of a friction term, thus 
providing an extra parameter y(a). It s also possible to generalize the model (1 1 3I > to mixtures of IMP ffl . 
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